Discrimination of the Healthy and Sick Cardiac Autonomic Nervous System by a New 

Wavelet Analysis of Heartbeat Intervals 
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We demonstrate that it is possible to distinguish with a complete certainty between healthy subjects 
and patients with various dysfunctions of the cardiac nervous system by way of multiresolutional 
wavelet transform of RR intervals. We repeated the study of Thurner et al on different ensemble of 
subjects. We show that reconstructed series using a filter which discards wavelet coefficients related 
with higher scales enables one to classify individuals for which the method otherwise is inconclusive. 
We suggest a delimiting diagnostic value of the standard deviation of the filtered, reconstructed RR 
interval time series in the range of ~ 0.035 (for the above mentioned filter), below which individuals 
are at risk. 
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I. INTRODUCTION 

Measurement of heart rate (HR) and evaluation of its 
rhythmicity have been used for a long time as a simple 
clinical indicator jjj., The main adaptive regulation of 
the sinus node function and thereby the HR, is exerted 
by the autonomic nervous system. The sinus node of 
the heart is a major organ in the integrated control of 
cardiovascular function. HR abnormality may therefore 
be an early or principle sign of disease or malfunction. 

Research from the last decade indicates that a quan- 
tification of the discrete beat to beat variations in HR - 
heart rate variability (HRV) may be used more directly to 
estimate efferent autonomic activity to the heart and the 
integrity of this cardiovascular control system Q). The 
finding that power spectral analysis of HRV could be used 
as a marker of cardiac autonomic outflow to the heart, 
was considered a breakthrough for clinical research ||Q . 

Autonomic dysfunction is an important factor in a 
number of conditions. In diabetes, an abnormality in 
autonomic nervous function signals an adverse progno- 
sis and risk of subsequent heart disease. Recognition of 
early dysfunction is therefore important. In overt heart 
disease autonomic imbalance is of significant importance 
in the pathophysiology of sudden cardiac death. Abnor- 
mal autonomic balance is an important prognostic factor. 
In heart failure this control system may be significantly 
deranged. 

Techniques which can discriminate the healthy HRV 
profile from a sick one are therefore highly desirable. So 
far this has not been accomplished, as a considerable 
overlap between healthy and sick, (i.e. healthy and dia- 
betes) H or high and low risk heart disease patients || , 
have been reported. The time series used for HRV analy- 
sis are derived from 24-hour ECG recordings. These are 



clinically widely used and offer important additional in- 
formation. However, several problems have limited the 
use and interpretation of the spectral analysis results. 
The ambulatory time segments inherently lack station- 
arity. Furthermore, they often include transients caused 
by artifacts, ectopic beats, noise, tape speed errors which 
may have significant impact on the power spectrum [0. 
This significantly limits the sensitivity of this technique, 
and thus may limit its applicability. 



II. METHODS 

One of the most successful techniques to analyze non 
stationary time series is the Multiresolution Wavelet 
Analysis This technique was recently utilized in 

order to analyze a sequence of RR intervals |l3|,[l4| . Ref. 
p3| identifies different scaling properties in healthy and 
sleep apnea patients. In a previous study, Peng et al Jig] 
were able to distinguish between healthy subjects and 
patients with heart failure by the use of the detrended 
fiactuation analysis. Later, Thurner et al jlj] used a 
similar procedure but focused on the values of the vari- 
ance rather than on the scaling exponent. For the scale 
windows of m = 4 and m = 5 heartbeats, the standard 
deviations of the wavelet coefficients for normal individ- 
uals and heart failure patients were divided into two dis- 
joint sets. In this way the authors of ref. Q succeeded 
to classify subjects from a test group as either belonging 
to the heart failure or the normal group, and that with 
a 100% accuracy. 

The Discrete Wavelet transform is a mathematical 
recipe acting on a data vector of length 2 m , m = 1, 2, . . . 
and transforming it into a different vector of the same 
length. It is based on recursive sums and differences of 
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the vector components; the sums can be compared with 
the low frequency amplitudes in the Fourier transform, 
and the differences with the high frequency amplitudes. 
It is similar to the Fourier transform in respect of or- 
thogonality and invertibility. The wavelets are the unit 
vectors i.e., they correspond to the sine and cosine ba- 
sis functions of the Fourier transform. One of the basic 
advantages of wavelets is that an event can be simultane- 
ously described in the frequency domain as well as in the 
time domain, unlike the usual Fourier transform where 
an event is accurately described either in the frequency 
or in the time domain. This difference allows a multi 
resolution analysis of data with different behaviour on 
different scales. This dual localization renders functions 
with intrinsic inaccuracies into reliable data when they 
are transformed into the wavelet domain. Large classes 
of biological data (such as ECG series and RR intervals) 
may be analysed by this method. 

Heart failure patients generally have very low HRV val- 
ues. To further explore the potential possibilities of the 
Multiresolutional Wavelet Analysis we have investigated 
a test group of 33 persons, 12 patients and 21 healthy 
subjects. The patient group consisted of 10 diabetic pa- 
tients which are otherwise healthy and without symp- 
toms or signs of heart disease, one patient which have 
had a myocardial infarction and one heart transplanted 
patient in whom the autonomic nerves to the heart have 
been cut. 
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FIG. 1. RR interval vs. (heart)beat number for a healthy 
subject. 

We have in the present study applied the same tech- 
nique as used in ref. [ [lit an d have by Multiresolution 
Wavelet Analysis been able to identify correctly all but 
one of 33 test persons as belonging to the group of healthy 
subjects or subjects suffering from myocardial infarction. 
The heart transplanted patient was included as a subject 
displaying the ultimative cardiac autonomic dysfunction 
- complete denervation. 



We have, however, elaborated on the procedure ap- 
plied in ref. fl4|| by utilizing a filter-technique. Thus we 
perform an Inverse Wavelet Transform, but retain only 
a specific scale in the reconstruction of the time series; 
a complete separation is observed for m = 4 or m = 5. 
In this way a reconstructed and filtered time series is 
obtained and a comparison with the original time se- 
ries shows a substantial difference in amplitude between 
sick/healthy subjects relative to the difference found in 
the original RR interval time series. The choice of m = 4 
or m = 5 was motivated by the findings in ref. Jl4| and 
by our own results. 

III. RESULTS 

We have calculated the standard deviation a wave for 
Daubechies 10-tap wavelet versus the scale m, 1 < m < 
10, for 33 persons. In accordance with ref. [Q we find 
that for 4 < m < 6 the a wave separate the two classes of 
subjects and hence provide a clinically significant mea- 
sure of the presence of cardiac autonomic dysfunction 
with a 97% sensitivity. This supports in a convincing 
way the findings of ref. [[t4| . We have been able to con- 
firm this trend with other wavelets. 

The main result of this study is however the possibil- 
ity to display the standard deviation of the RR inter- 
val amplitude vs. the beat number in the reconstructed, 
filtered time series. This standard deviation, here de- 
noted by a filters can be used to obtain a separation of 
sick/healthy subjects. 

In fig. [I] we display the RR intervals vs. the beat num- 
ber of a normal subject. The wavelet technique cleans 
the highest and lowest frequencies from the overall pic- 
ture. The highest frequencies contain noise and the low- 
est frequencies contain mainly external influences on the 
HR pattern like movement and slower trends in HR level, 
which are not necessarily reflective of autonomic nervous 
activity. After the removal of these frequencies one is left 
with the characteristic frequencies of the heart. 

Fig. H shows the standard deviation a wave for a 
Daubechies 10-tap wavelet as a function of the scale num- 
ber m. The almost total separation between sick and 
healthy subjects is obvious. 

Patient #1, falling into the range of sick patients, has 
a very low HRV both on a 24-hour scale and short term. 
The patient is a survivor of a heart infarct and is at high 
risk of sudden cardiac death. 

Patient #2 has the lowest a wave values in the range 
4 < m < 6. He has undergone a heart transplant; the 
nerves to the heart have been disconnected and there is 
almost no HRV. 

Patient #3 is a diabetic patient, who is classified by 
the wavelet technique as a high risk patient. Diabetic 
patients with abnormal cardiac autonomic function have 
an adverse prognosis and increased risk of heart disease. 

Patient #4, also a diabetic, seems to be less at risk. 
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His 

Pwave is near the transition between healthy and sick 
subjects. 
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FIG. 2. Daubechies 10-tap wavelet. a waV e, the standard 
deviation, is plotted as a function of the scale m, 1 < m < 10. 
The corresponding window size is 2 m . The empty symbols 
indicate the healthy subjects, the opaque symbols indicate 
patients. The circles designate normal subjects, the squares - 
diabetic patients, diamond - patient at risk with heart infarct 
and triangle - a heart transplanted patient. 



individual hints to the possibility that the correspond- 
ing frequencies are characteristic of those at which the 
autonomic nervous system works. 
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The method used in ref. |l<|] fails for subject #5, who 
appears in the risk group, although he had no evidence 
of diabetes or heart disease. 
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FIG. 3. Daubechies 10-tap wavelet filtered inverse trans- 
form. The symbols are as in fig. 2. 

In fig. H the standard deviation of the amplitude 
of the reconstructed time series has been calculated for 
1 < m < 10. Again, a total separation between sick and 
healthy subjects is apparent. The fact that the (Tfuter re- 
main almost constant for scales between 4 and 6 for each 
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FIG. 4. (a) Typical time series segments for a sick and 
a normal individual, (b) Typical reconstructed, filtered time 
series for the above individuals. The segments shown are the 
same as in (a). The filter is created by the inverse transform 
of coefficients with scale m — 4. 

Fig. |^a shows a typical RR interval time series for a 
healthy and a sick subject, whereas fig. ||b shows the 
reconstructed time series (m = 4). One notices that the 
difference in amplitudes for healthy/sick subjects is much 
more pronounced in the latter time series. 



3 




2000 3000 4000 
Fourier index 



6000 



risk. 



0.10 



is now properly classified as not being at 



0.08 



P 0.06 
i 



0.04 



0.02 



0.00 



O healthy 

□ diabetic 

O heart infarct 

A heart transplanted 




6000 



, Fourier index 

The Fourier transforms of the above 



FIG. 5. (a) and (b) 
(fig. ^). An index of 1000 represents a frequency of 0.02 Hz. 

Figs. ||a and ||b show the Fourier transforms for the 
time series displayed in figs. [|a and respectively. 
These power spectra appear similar, however differ in 
their respective order of magnitude. Clearly, the recon- 
structed filtered time series are distinct by the amplitude 
as well as the broadness of their Fourier transforms. 

In fig. H we have obtained a complete separation be- 
tween the sick and healthy subjects by application of a 
filter which is created by retaining wavelet coefficients 
with scales 1 < m < 6. This filter was motivated by the 
observation that a separation is evident for these scales 
(see figs. Hand||). One observes that the healthy subject 
#5, who failed the wavelet transform diagnostics of ref. 



FIG. 6. Daubechies 10-tap wavelet filtered inverse trans- 
form. The symbols are as in figs ^| The filter is created by 
the inverse transform of coefficients with 1 < m < 7. 



IV. CONCLUSION 

Our study supports the conjecture of ref. H] that 
healthy subjects exhibit greater fluctuations (larger 
&wave values) than patients. This difference in fluctu- 
ations become most evident on the scale 4 to 5 (corre- 
sponding to windows of 16 and 32 heartbeats), but in our 
study it is apparent at all scales from 1 to 7 (windows of 
2 to 128 heartbeats). 

The most distinct difference between sick and healthy 
individuals appears in the amplitude changes in the 're- 
constructed' time series, where the windows of 16, 32 
and 64 heartbeats contribute in a similar way. Letting 
the window be as small as 2 4 heartbeats is enough to 
allow the healthy group to show substantial variation in 
the size of RR intervals implying a large a value, but 
is at the same time too small a window to let the sick 
cardiac autonomic nervous system introduce significant 
variations in the length of the RR intervals and hence 
allows it only to reach a a value essentially smaller than 
the healthy heart. 

The final conclusion of this study is that in order to 
obtain a complete separation between healthy subjects 
and patients one has to consider a range of scales (as 
shown in fig. ||) instead of only one scale (as in figs. |2| 
and ||). This implies that, a filter as m fig- H can be 
used as a diagnostic indicator, with a delimiting value of 
~ 0.035 (for the above mentioned filter), below which the 
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persons have abnormal cardiac autonomic function and 
will be at risk. 
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